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APSTRACT 

Geostropnic relative vorticity is commuted along three latitide 
circles reneralli-:y encompassing the United States, using a Fourier analy- 
sis of the 500=-mb neicht field for the veriod 3-6 November 1961. 

“ statistical comynarison of vorticity obtained in this manner with 
it*s socuare grid computation is mede. 

The writers are deeply indebted to Professor F. L. Martin of the 
United States Neval Postgraduate School for his suggestion of the topic 
and his continued kelp throughout the investigation and during the 


preparation of this paper. 
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lene roane 10, 

The advent of the electronic computer coupled with increased data 
coverare hag removed the primary obstacles to the use of harmonic analy- 
sis as a useful instrument in atmosrheric research. 

Use is made of a Fourier analysis of the 500-mb height field. A 
computer program devised by G. Arnason and i, Reese.and made available 
by Fleet Numerical Weather Facility! was used to obtain the harmonic 
analysis sround latitude circles at distances of 1 through 27 grid lengths 
d = 381 km from the pole. An equation for relative vorticity using the 
results of the spherical harmonic analysis is derived and tested statisti- 
cally against the more usual numerical computation of reostrophic relative 
vorticity made using a square grid. 

The problem undertaken was to compute relative vorticity at 500 mb 


using the geostrophic vorticity equation in its spherical coordinate form 
Ha, p. 358. 
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In equation he), Z was replaced by its expression in terms of a Fourier 
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Here z is the zonally averaged contour height, whereas z* is the perturba- 
tion quantity which remains when the z field is subtrected from the actual 
height field. The computer printovt included maps of each of these three 
fields at 500 mb for the period under investigation, 3 through 6 November 
1961. Geostrophic relative vorticity was computed by the above mentioned 


technique for 17 longitudes along each of three latitude circles generally 


1 ; | 
Henceforth referred to as FNWEF 








encompassing the north-south extent of the United States. This gave a 
total samnle size of 204 computed relative vorticities which were then 
correlated with square-grid computations of relative vorticities obtained 
fron FNWF charts of absolute vorticity for exactly the same set cf points 
and times. 

2. The geostrophic relative vorticity in terms of a harmonically analyzed 
contour field. 

The adequacy of the Fourier representation of meteorological pare- 
meters will not be discussed here; the interested reader is referred to 
aie 

It was previously mentioned that use was made of the FNWF computer 
orogram for a Fourier analysis of the 500 mb contour field. A brief de- 
scription of the theory and mmerical evaluation of the Fourier coeffi- 
cients is given in (3) and only the applicable portions of the program 
will be mentioned here. 


The zonal mean height 2 is defined as 
_ a 
mf - ¢& 
Pais ane &D Ey) 
and its zonal perturbation, in accordance with equaticn (2), is 
ee oy Me 
ee (4) 


Expanding 2 in a Fourier series one obtains at a given latitude 
Z eee) AP Vee Ses ea i (5) 


where A and B are the coefficients of the nth harmonic functions in the - 
Fourier expansion of 2, and Mis the total number of harmonic waves neces- 


sary to resolve the contour field. Replacing A, and B, in (5) by an am- 
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plitude factor C. and a phase angle = , (5) may be rewritten as 
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Equation (6) has been vrogrammed for the CDC=1604 dipital computer by 
FNWF, and the program includes: a contour grid map of the zonal pertur~ 
bation field z*, and a tabular list of the Fourier parameters C,, ©, and 
2 for each of 27 latitude circles. The last three paremeters are actually 
computed in relation to a zonal band of width equal to the grid mesh d of 
the FNWEF map grid, but may be considered applicable at the center of the 
band, 

Using soherical coordinates \ and ¢ (longitude and latitude, respec~ 


tively) the derivative operators 
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If now z is replaced in (7) by its exnansion in a Fourier series as given 


in equation (2) we obtain 
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Remembering that C and * are independent of longitude )) but not of 





latitude }, we may perform the indicated differentiations and combine 


terms, with the result that J, in spherical coordinates becomes 
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Equation (9} may be shown to have the form 
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Were the superscript bar denotes an average with respect to n over all 
eapmificant wavés n= 1, ....., M. 

The FNWF maps are stereographic maps true at latitude 60N, and else- 
where the grid mesh, d = 381 km, is related to the true distance d_ by the 
relationship d = md_, where m = (1 + sin 60)/(1 + sin )) and is the so- 
called map factor, 

The partial derivatives in (11) may be evaluated by centered finite 


differences, The numerical values of C, and %, were simply taken from 
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the tabnlar list of computer output of the harmonic analysis at the vari- 
ous latitude circles. For example, 3-, /; may be written using 
finite differences between adjacent latitude circles (for which aa¢  d/m) 


in the form 
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The latitude circles of the Fourier analysis employed in the computa- 
aon of Sg corresponded to N = 12, 15, and 18 grid mesh distances from the 
pole, which in turn corresnorded to the latitude circles: 4A8.9N, 38.6N, 
and 28,3N (see fig. 1). However, in the computation of first and second 
derivatives of the nth Fourier wave parameters, these paremeters at lati- 
tude circles N=-71 and +1 are also needed. 

The term d2/ A is more tractible in the present case if the finite- 
difference substitution r{\ = 24/n is made; thus, the centered deriva- 
tive dH AN in finite differences, assumes the form 
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where “= ¢cG-: $ is the radius of the latitude circle, and the subscript 
mri is one grid mesh eastward of the ith point along the latitude circles 
Including the form of !“"°/3A given in (12), equation (11) may now 


be written in the more compact form 
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The averaging with respect to n in the expressions for J) and W, is 


| 
assumed to be valid, accordine to the mean value theorem. Actually how- 
ever the average was approximated by a simple arithmetic mean of the con- 
eribievons [or n= 13ige.., M. Inus it Waemnecegsary to canpnve contri— 
butions for each n and divide the result by the number of waves, 

ean Mieghem [ 2 | noints out that waves of space number higher than 12 
to 17 contribute little to the variance of the contour-height field. Con- 
sequently, the Fourier analysis was truncated at M = 12 at the two higher 
latitide circles, | although wave numbers 1 through 15 were employed at the 
lovest latitude. 


The coefficients ‘J, and W, are constants for a particular time and 


latitude circle, and can be computed directly from the tabulated harmonic 
data, These values are listed in table 1. However, in applying equation 
(13) it was necessary to internolate values of z2* from the perturbation 
map at the points indicated in fig. 1 along the specified latitude circles. 
These roints were arbitrarily chosen at integral multiples of five degrees 
of longitude, encompassing the west-east extent of the United States. 
Therefore, for ease in selecting data points, a grid overlay was con- 
structed presenting all latitude circles N=11,...-,19, drawn so as to in- 


corporate the map factor, 


Tat the two higher latitudes, ‘!, was negative ensuring thet lows 
corresponded to positive relative vorticity. At the lowest latitude 
(1=18), it was necessary to use M=15 in order to ensure that this condi- 
tion was met, 





TABLE 1: Values of W, ond WV, to be applied in equation (; >. Al] values 





me 
are in Wnagewo: 10 See My - 
eS ee Ss LATITUDE CIRCIE _ aa 
ena os ee Se ae eee | pea 
DATE N= 12 | N=15 —  N=18 
a Sr y= -.1467 Woo= —.1635 Woo= -,2696 
3 Nov. 1961 
| W =-,0079 | W = -,0370 Wo = =.0303 
| or eee ee eee ee 
| Wo =-.1490 | UW = 0617 W = -.0420 
4, Nov. 1961 | | 
Wo= +0311 § W = ,9316 W = 0138 
—— mh = 13556 WW) = 11012 i See TES 
5 Nov. 1961 
Wo = -,9308 Wo=-.0230 i W = ,1682 
—_— “oY 2 oss”  Welwato vel eae Oo oe 
= © Nov 64061 
W = 0757 W = 0349 W = -,761 O 
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Note that the term (ZU, *: 4¢)/ in both (11) and (13) also varies 
from point to point around any latitude circle. However, the variation of 
this term was so small that the mean zonal wind for the latitude was used. 
The mean zonal wind Us was computed using the latitudinal gradient of 


zonal mean height z, that is from 


again using finite differences, 

Along each of the three latitude circles, 17 values of W 2* and 
W iy 2%, ) were obtained for each of the four days, giving a total of 
204, separate values of these quantities. 
peweotatistical results. 

For precisely the same set of voints and davs as described above, 
corresponding values of 5. were obtained from grid maps of geostrophic 
relative vorticity made available by FNWF, Their maps essentially pre- 


sented the field of Jy computed using their square-grid mesh according 
8 





to the usual fornula 


! up S 
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where Z is the four-point space-mean contour height and ad = 381 km at 
latitude 60N, 

Since the averaging process in equation (11) is somewhat approximate, 
it was felt that better verification would result if the relationship be- 


tween | sand Wat and W. (2# -2% ) was written 
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Here &,, A, and 4, are regression coefficients" to be determined by a 
least squares technique applied to the 204 simultaneous values of Wize 


and W. (2% ~2*, ). Because of occasional large gradients ones over 5° 


WF 


longitude intervals, the values oa eevee 


were smoothed by the three-point 


binomial running mean formula 
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where X_ is any value of a interpolated from the FN vorticity-grio 
map at each point. The regression equation actually involves the correla- 


tion of the smoothed } values with the independent variables already 


ad 


noted. 

Using the BIMD 06 statistical program as adapted to the CDC~1604 
computer, it was necessary to type 204 IBM cards containing the smoothed 
J wwe Values and the simultaneous values of the independent variables of 
the right side of equation (15), Best-fitting values of the regression 
coefficients were 


Itrese regression coefficients are non-dimensional since the units 
of W 2° and ez.) are in sec™!. 





Ay on = .1719, A, = OSS, and As = 0.50524, 


1 and the multinzle linear correlation co- 


The standard error of estimate 
efficient for the regression equation were 20.3523 and 0.7632 respectively. 

Panofsky and Brier | 4, ile 113] present an F test for the analysis of 
the variance explained by the resression equation, The F value may be de- 
fined by the ratio of the explained variance to the residual variance 
relative to the regression equation. This formula is 

Fe RPh) wae 
CLs Ie ) ne 

With p = 2, the number of independent variables, and’) = 204, the number 
of random samples, the value of Ff from equation (16) turns out to be 
139.33. Nowever, assuming normally distributed independent and dependent 
variables, the critical F at the 9.99 significance level. with 2 and 201 
decrees of freedom is 4.71. Thus the regression equation (15) is signi- 
ficant at a very high confidence level. 


The partial correlation coefficients of W 2* and Witcgaten ) were 


= ae 
0.7620 and 0.0382, respectively, indicating very little contribution to 
the regression equation for the latter variable. An intuitive explana- 
tion for the insignificant contribution of Wy (2% —2% ) may be obtained 


from fig. 2. The differences, 2' “Zi taken over the pressure centers 


would be zero, whereas in actuality the largest absolute values of Jy 


would be found at tnese points. Ilikewise for a sinusoidal contour wave 


ae 
ray 


pattern maximun magnitudes of 2* -2%* wonld exist near contour inflection 


A045 A! 
points, where J4 .~°- 
The possibility of "automatic" correlation, viz, correlation between 


the two independent variables, was examined and none was found to exist, 


1 . = = 
All values of ‘xy are in units x 10 > sec ue 
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The individual latitude circles vere statistically examined in a like 
manner. The two northernmost latitude circles gave a multiple linear cor- 
relation coefficient of 0.84, whereas there is no significant correlation 
between See and the two independent variables W,2?* and Wa (Zs =2%) at 
latitude 28.3. This inconsistency at low latitudes is attributable to 
"noise" in the Yourier renresentation of the contour field. The "noise" 
is due to the smal] variance of the height field at this latitude, coupled 
with the relatively large RMS error possible in the 590 mb analysis in the 
subtropics during the veriod of study. The use of second derivatives of 
the Fourier parameters also increases the amount of "noise" generated in 

he determination of F : 
4. Conclusions. 

Theoretically trere s*onuld have heen a large positive multiple linear 
correlation between Vs and Wo2* and a A cee This proved to be the 
case at latitudes 4&.9N and 35.6N where the contour field was well resolved 
by twelve harmonic waves. However, at the lowest latitude circle, 28.3N, 
the several sources of errors listed at the end of section 3 could, and 


probably did, influence the low correlations that were obtained. These 


were basically due to the lack of resolution of the contour field by the 


We 





Fourier analysis. However, t»e same reasons would cause the square-srid 
commitation of J to be subiect to large errors. 


Based upon the results at the two higher latitudes, it seems a 


4 


pparent 
that the use of a Fourier analysis of a meteorological field presents a 
very valuable tool for atmospheric research, With the advent of the 
electronic computer and machine-analyzed contour charts, the use of har- 
monic analysis is an effective and economically feasible avproach to the 


computation of synoptic aids for the practising meteorologist. 
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